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16. Abttroct 

The problem of calculating the Inlrlal distribution of mag- 
netic induction in the armature stationary relative to the p 

system of a hysteresis clutch is considered. Using several 
assumptions, the problem Is reduced to calculating the statl 
magnetic field In the ferromagnetic plate with finite and 
continuous magnetic permeability placed In the air gap betwe 
two Identical, parallel semiconductors with rack fixed rela- 
tive to the "tooth" or "slot" position. 13 ref. In Russian. 
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CALCULATION OP THE MAGNETIC FIELD IN THE ACTIVE ZONE OP A 

HYSTERISIS CLUTCH 

M. A. Ermilov, Candidate, Technical Sciences 
0. M. Glukhov, Engineer, Moscow Aviation Institute 


Hysterlsls clutches are a widely distributed type of asynchronous- /6l* 
synchronous clutch with magnetic coupling. Structurally, a hysterlsls 
clutch is a dual rotation device. On one of the rotating parts - usu- 
ally the driving part - there is a pole system of the Inductor. On the 
other rotating part there is an armature in the form of a cylindrical 
or disk-type layer of ‘'hysterlsls” material (l.e., solid magentlc ma- 
terial with higher- than- usual remagnetlzatlon losses). Transmission of 
torque from the hysterlsls clutch is the result of the interaction be- 
tween the magnetic field of the Inductor pole system and the field due 
to the residual magnetism of the hysterlsls layer. If the hysterlsls 
layer consists of electrically insulated rings made of solid magnetic 
material, then the transmission of torque from the clutch is practically 
Independent of slippage [1-3]. Hysterlsls clutches (HC) which have 
permanent magnets are called magneto-hysterlsls clutches CMC) in con- 
trast to clutches with electromagnetic excitation (EHC). 

The static and dynamic characteristics of EHC are analogous to the 
corresponding characteristics of electromagnetic particle clutches (EPC) . 
With respect to over-all mass and energy indices, the ECH are inferior 
to the EPC; however, they have higher reliability indices and stable 
characteristics. The ability of the EHC to transmit torque proportional 
to the excitation current permits its use in high-speed servomechanisms 
with an actuating device operating on the controlled electromagnetic 
clutches of a transmission [4], Prom the point of view of high-speed 
action and control capacity, the most efficient design for an EHC in- 
cludes a driven rotor in the form of a thin walled sleeve and a "bomb- 
type” system of poles with a two-way, "checkerboard" distribution of 

•Numbers in margin Indicate pagination in foreign text. 
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the toothed poles (relative to the walls of the sleeve armature) [5], 
Flfrure 1 contains a schematic drawlnp; of an EHC of this type with a 
stationary field wlndinfr [6]. ^ Is the drive shaft; 2 Is the driven 

shaft, 2 Is the stationary mapinetlc circuit; ^ the excitation coll, 

5 Is the external magnetic cl.'cult, 6 Is the external polo ring, 7 
Is the .junction element made of non ferromagnetic material, 8 Is the 
driven rotor (armature), 9 is the Internal pole ring, the dotted line 
shows the path of the basic magnetic flux. 



Figure 1 . 

The procedure for calculating the maximum torque delivered by an EHC 
for a given control current when the drive shaft and the slave shaft 
are rotating In synchrony Is based on the so-called energy approach to 
the analysis of the operation of hysterlsls electrical machines [2J. 

In accordance with this approach, there Is an "Initial" distribution 
of magnetic Induction Into the hysterlsls layer, l.e., after the field 
winding Is first plugged Into the constant voltage and the ’urrent Is 
established In the winding, then the stationary magnetic field Is dis- 
tributed In the armature, which Is stationary relative to the polar 
system, and applied to the EHC. Linder the conditions stated, the mag- 
netic state of the armature Is characterized by the points on the Ini- 
tial magnetization curve of the solid magnetic material employed. On 
the basis of the Initial Induction distribution, the parameters of the 
remagnetlzatlon cycle are found for each elementary ring layer of the 
solid magnetic material. These determinations are made when the polar 
system Is rotating at a very low speed relative to the braked armature 
- a speed which ..s insuf '’1 dent to cause eddy currents to appear in the 
armature. In conformity with the energy balance, the electromagnetic 
moment Is directly proportional to the energy losses during a single 
remagrietlzatlon of all elements of the armature ring layers. 
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An approximate solution of the problem concernlnR the Initial dis- 
tribution of Induction In an EHC of the typo under consideration Is pre- 
sented In [7] and [8], In this case, a technique for simplifying cer- 
tain boundary value problems In magnetostatics Is used which Is analo- 
gous to that In [9]. However, such Important questions as the' In- 
fluence of the extent of overlap of the toothed poles, and the effect 
of the slzs of the air gap between the armature and the system of poles 
has on the field pattern ir the armature have not been Investigated. 

In the present article, we present a more rigorous solution of the 
problem concerning the Initial distribution of Induction In the active 
zone of an EHC. This makes It possible. In particular, to evaluate 
the effect produced by the degree of overlap of the toothed poles and 
by the size of the operating air gap on the diagrams of the vector 
components of the magnetic induction in the armature. 

This solution was obtained under the following assumptions: 1) the 
depth of the operating gap between the toothed surfaces of the polar 
rings Is small Incomparlson to the length of the active zone In the 
direction of the rotation axis of the clutch; therefore the field In the 
active zone Is plane-parallel; 2) the toothed spaclngs of the polar 
rings are small compared with the diameters of the operating annular 
gap; hence, the cylindrical surfaces of the armature and the crowns of 
the teeth can be replaced by parallel planes (of Infinite extent); 

3) the teeth along both sides of the armature are rectangular In shape, 
and have Identical dimensions and spaclngs, 4) In the elements of the 
active zone there are no eddy or displacement currents (since the state 
of the magnetic system Is under consideration, when the armature Is 
stationary relative to the toothed poles and the excitation current Is 
constant after the winding Is first plugged Into the constant voltage); 

5) the magnetic permeability of the solid magnetic armature material 
does not depend on the field Intensity and has a finite value which Is 
determined from a linear approximation of the magnetization curve; 

6) the n. 'gnetically soft material of the teeth Is unsaturated and thus 
its magnetic permeability Is Infinite. 

With the aid of the above assumptions, the problem Is reduced to 
calculating the magnetic field In an Idealized model of the active 
zone; l.e., calculating the magnetic field in a magnetic system of two 
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Identical parallel semiconductors having the profile of a toothed rack 
and fixed in a ’’tooth opposite gap” relative position; in the gap be- 
tween these two semiconductors there is placed a plate of ferromagne- 
tic material with constant, finite magnetic permeability. 


Figure 2 shows a partial cross-section of the model of the active 
zone. The width is equal to that of the toothed spacing. The follow- 
ing notation is used for the geometrical measurements: denotes the 

toothed spacing, a, the width of the slot, h, the depth of the slot 
(height of a tooth); the thickness of the ferromagnetic plate, 
representing the armature; _6, the size of the unidirectional gap. 


The plane-parallel magnetic field in the region under consideration /62 
satisfies the Laplace equation for the scaler magnetic potential, 
u(x, y), when no current flows in the region: 




( 1 ) 


The equations 

da 

du> 

liy — — fjy » 


( 2 ) 


express the relationships between the projections of the magnetic in- 
duction vector and the potential function, is the magnetic con- 

stant, is the relative magnetic permeability of the region. 



Figure 2. 

A) Region 3 B) Region 2 C) Region 1 D) Region 4) E) Region 5 
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Analysis of the probable paths of the magnetic flux in the model of 


the active zone (PlRure 2) Indicates the presence of periodicity and 
symmetry properties in the potential function. The period in the longi- 
tudinal direction equals the toothed spacing and is identical for the 
upper and lower semiconductor, hence 


6 , 


/ A \ 


(3) 


The pattern of the lines of force of the field under study is char- 
acterized by plane symmetry with respect to the trace of the plane 
X ■ 0, and by central symmetry relative to the points A(t 0) and 

Q(-t /^, 0) . Plane symmetry is due to the fact that the semlcon- 

z 

doctors are fixed in the "tooth opposite gap" relative position, and 
central symmetry results from the fact that the profiles of the toothed 
surfaces of the semiconductors are identical. Mathematically the pro- 
perty of plane symmetry is expressed by 

/«, (X. //) ^ — n, (-*, i>): I 


or 


H (X, y) •» (I (—X. y). 

The central symmetry property is expressed by 

II (A. y) » It A. — 


or 


tt(x. 1 ^) « • m^4/— X, -yjforu^~l~i oj-o. 


(5) 


( 6 ) 


From the plane symmetry and periodicity properties. It 
the axes of the teeth (slots) coincide with the magnetic 


lines. Therefore, ,tn i a / a . 

'’“I A ^ / A , . , \ 

‘'for v*' (V+^ ‘ ''j' 


I'ollows that 
induction 


(7) 


and the least recurrent part of the active zone is Included between 

the lines x * 0 and x ■ t /2. This part consists of regions with dlf- 

z 

ferent magnetic permeability, l.e., it is piecewise uniform. The ana- 
lytic expression for the potential function In each of the regions can 
be found by realizing one of the methods for solving boundary value 
problems for the Laplace equation. In [10], an effective method Is 
reconjnended for finding the potential function In an air region with 
orthogonal boundaries when the potential varies on the boundarv of the 
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atr gap and the ferromagnetic plate according to a known law. This 
method Is based on representing the potential function In a trigono- 
metric series In the slot and the air gap, and subsequently connecting 
the solutions at the separation boundary between the slot and the gap 
by means of an Infinite system of equations. This method, which Is 
conventionally designated as "sewing together", has been successfully 
applied by a number of Investigators to the calculation of the magnetic 
field In the active zone of Induction type electrical machines. Below, 
we present Its realization as applied to the formulation of the above 
problem . 


The region which models the EHC active zone Is divided Into five 
rectangular subregions with sides parallel to the coordinate axes: 
the ferromagnetic sheet, the air gaps 2 and the slots 3 and 5 
(Figure 2). The external boundary conditions for calculating the field 
In this system of subregions are the conditions (7) and the equl-po- 
tentlallty condition for the toothed surface of a semiconductor. The 
difference In the potential of the semiconductors Is equal to the mag- 
netizing force, expanded In conducting the magnetic flux across 

the operating air gaps and the ferromagnetic plate. 

From central symmetry (6) It follows that when calculating the po- 
tentials relative to the point A (l.e., when u^ * 0), the potentials 
of the upper and lower semiconductors are equal modulo and 

are opposite In sign. 


For the purpose of representing the results of the solution In cri- 
teria form, vie Introduce relative values of the coordinates, the dimen- 
sions, the potential and the Induction; as basic values for the length, 
the potential and the Induction we adopt 

^ Q,oF ^ 

/«•=>/.; ■7^ ** (S J • 


The relative values are connected by the pertinent relationships 

• i 

/(lUrnx, y, a, h. A) -~7~> 




Hereafter, the sign * will be omitted. 


6 


In accordance with the recommendations In [11], the potential func- 
tion on the boundary of regions 2 and 2 1® represented in the form 

Ob 

( 2/1 — 1 ) 

Z (8) 

MM-i. i. .. 

_ d ll 

for — 'y < A 

for the realization of the "sewing together" method. 

The trigonometric series In (8) has no even harmonics; however, It 
was decided not to take account of this fact In the numbering of the 
amplitudes of the harmonics, c^^, since by so doing, the programming of 
the problem for a small ilgltal computer Is simplified. In particular, 
the formula for computing the coefficients of the system of linear al- 
gebraic equations obtained by "sewing together" subregions are described 
more compactly. Tne matrix representation of this system of equations 
Is 

(M. ^9) 

where m Is the number of rows; n the number of columns, ^], the 

square coefficient matrix of the system, [c^^] Is the column matrix of 
unknowns]; CSjjj] Is the column matrix of the right-hand members of the 
system. 

The scheme for realization of the "sewing together" method Is de- 
scribed in sufficient detail in [12]. With the aid of an analogous 
scheme, the following formulas for computing the elements of the matrix 
were obtained: 


, ttrJjZ . ( 10 ) 

II (IM I {>m — if* 

n» f- (Jm — I) {'in — I) 

00 

n. i(-i ) • • I f h 'ikr.i I t- »h ^ ^ ^ ^ 

K-h* 1 chV/lwJJ th -ikr.ii + »h rinA * 


where 
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Ym« m ^ 


when m « n 
when m K n 


( 12 ) 


ft,*"- 


(—1)* COS ktui _ 2/ — I 




when 


/63 


(13) 


21- I 
2k ‘ 


’It, 


Also, in (13) 1 identifies the first subscript - i.e., during the cal- 
culation, substitute m or n for i. 


For the potential function in the region 1, the following expression 
is obtained; 


«(<. — oISa+Th? 


* + « 2j 2n.— 1 + 

n=:i, 2.3,,.. J 


00 00 


/».*fsh2to*(i, + 0,5A) + 
JL^Wr -K-l)*»h2fe»(»-0.SAH . ^ 

Zilj f(— O'* + Ch 2*«AJ sh 2kr.8 + 

* ‘ ' +*h2fctA ch2teS 


(1^) 


AbbI A:bs| 


On the basis of (2) and (l4), the projections of the magnetic induc- 
tion vector in region 1 are determined:. 


, " ■ . 00' 00 

, 8h 5kr. (!/ 4- 0,SA) + (—!)>» sh 2Ax (y — 0>5A) „ 

^ I(— 0* + ch 2 AbA) sh 2An3 -f- l*r' 


00 00 

+ 2«a2J2Jc,f,.KftX 

> ch 2Hm (tf + 0.5A) + (-1)> ch 2tm (y - 0.5A) 

^ I(— O’* + ch 2/nrA] sh 2 Ar3 + 1*-' sh 2AttA ch 2kr.8 ‘ 


(15) 


( 16 ) 


The average value (constant component) of the component of the in- 
duction vector normal to the surface of the ferromagnetic plate is the 
same in all horizontal sections of the regions 2, and 4: 
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I 


( 17 ) 


fO.S f m I 

ft frir • 1 L 

-“O.J 


This quantity Is connected to the magnetic conductivity of the 
flux In a toothed division (on a unit length In the active zone, direc- 
ted toward the axis of rotation of the clutch) by an equation which Is 
analogous to a familiar equation In the theory of Induction-type elec- 
trical machines [13]: 


where 


. u / 

» -I — I wb/a* OH 


• . o.v,* . 

•• 


(IP) 


(19) 


Also In (l8) and (19), and In all succeeding equations we use absolute 
units of measurement. Prom (17)-(19), follows: 


where 



( 20 ) 


( 21 ) 


Is the magnetic conductivity when there are no slots (l.e., v'hen 
the semiconductors are smooth): 


Analysis of the formulas for the coefficients of the system (9) shows 
that undue Increase In the depth of the Inductor slots Is Inadvisable 
In practice. In fact. It follows from (12), that when h > a the mag- 
nitude of Y changes very little as h Increases (since cth (2m-l)ir 

m , n — 

» 1, even when m » 1 , l.e., the field pattern preserves its previous 
form. 


From (10), (11), (15), and (l6). It follows that when the conditio. is 
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th '2knl 

are satisfied simultaneously, the ferromagnetic plate (armature) Is 
equivalent to a plate with Infinite magnetic permeability, and the 
field pattern In the air gap between the armature surface and the toothed 
surface of the semiconductor Is analogous to the pattern of the exci- 
tation field In the operating gap of an Inductor generator with constant 
flux [13]* Here, the Initial field distribution In the mass of the 
plate (armature) Is practically Independent of the values of Its mag- 
netic permeability. When A ^ 6, It Is sufficient to set k ■ 1 In (2U). 

The relations obtained were used In the formulation of the algo- 
rithm for computing the Initial Induction distribution In the EHC arma- 
ture. On the basis of this algorithm, an operating program for the 

digital computer "Nalra-2" was drawn up and calculations of the para- 

« « « « 

meter's a, h. A, £, and were carried out. The results of the calcu- 
lations were checked on physical electromagnetic models - In particular, 
on a large-scale mockup of the active zone of Inductor generators with 
ripple flux [13]. 

The ferromagnetic plate which models the armature was made from 
mark EKHZ chrome steel. This steel Is Isotropic with respect to mag- 
netic propei'tles. The maximum value of the sine of the hysterlsls 
angle for this steel Is attained during remagnetlzatlon along the par- 
tial hysterlsls loop with vertex at the point H ■ 6 ka/m, B ■ 0.9 

^ my my 

tl [2]. 

The relative magnetic permeability matching this point of the mag- 
netization cu.”ve Is ■ 120 relative units, also when < 1.1 tl, a 
straight line through the origin is a good approxlmaclon to the mag- 
netization curve. 


In one of the experiments on the model, the following values for 



the geometric dimensions of the active zone v^lements were obtained: 
t *70 mm, a ■ 35 mm, h • 35 mm, A ■ mm, ^ ■ 2 rrm. To these 
values there corresponds the following set of Initial parameters for a 

« tf II 

calculation In relative units: a ■ 0.5, h • 0.5, A • 0.0571^, and 

A • 0.02857. 

Figure 3 shows the curves of the Initial distribution of the tan- 
gent to the sui'face of the hyst oriels layer for the component of the 
vector of the magnetic Induction In the layer mass when ^ ■ 0 (curve 1) 
and at Its surface when ^ ■ 0.5A (curve 2). Those curves were construc- 
ted from the Initial data Introduced above. During the <“xperiment, the 
magnetization force, was fixed at 1115 a; also the basic value of 

the induction corresponding to 1 relative unit amounted to 0.01 tl. 

With the help of mllllweber meters, the normal component of the Induc- 
tion vector In the air gap was measured at the surface of the hysterlsls 
layer. 

Figure ^ shows the curves of the Initial Impact distribution nor- 
mal to the surface of the hysterlsls layer for the component of the mag- 
magnetlc Induction vector In the layer mass when ^ * 0 (curve 1) and 
at Its surfaci (curve 2). The curves were constructed using a digital 
computer; however, values of the Induction which were measured are 
Indicated by small circles. The high degree of linearity of the mag- 
netization curve of the hysterlsls material In the operating Induction 
zone (up to the region of the bend) contributed to the closi agreement 
between the computed and experimental data 

The analytic relations presented In the article show that quite /6^ 
a large number of parameters Influence the Induction distribution In 
the active zone of a clutch. A detailed Investigation of the Influence 
of each parameter can be conducted by Introducing statistical methods 
In planning a numerical experiment. The results obtained can be used 
when designing EHC, and also as a first approximation when computing 
the field In the active zone, taking into account the saturation of 
the toothed poles and the nonlinearity of the magnetization curve for 
the hysterlsls material. 
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